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Abstract 

Several upper bounds on the size of quantum codes are derived us- 
ing the linear programming approach. These bounds are strengthened 
Ch ' for the linear quantum codes. 

y. '. 1 Introduction 

> . 

^ , Recently P.Shor presented a polynomial time algorithm for factoring large 

^ I numbers on a quantum computer [19|. After this the interest in quantum 

computations grew dramatically. One of the crucial problems in implemen- 
tation of quantum computer appeared to be the one of eliminating errors 
caused by decoherence and inaccuracy. Unlike the classical information, the 
quantum information can not be duplicated Since error correct- 

ing codes protect classical information by duplicating it, their application 
for quantum information protection seemed to be impossible. However, in 
]l8[ P.Shor has shown that quantum error correction codes do exist and 
presented the first example of such one-error correcting code encoding one 
qubit to nine qubits. In M.Knill and R. Laflamme formulated necessary 
and sufficient conditions for an error to be detectable by a given quantum 
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code, and thus introduced the notion of the minimum distance of a quantum 
code. In [^] P.Shor and R.Laflamme showed that, similarly to the classical 
codes, the quantum codes have enumerators related by the MacWilliams 
identities. Properties of quantum enumerators were extensively studied by 
E. Rains |14],|15],|1(:]. In particular, he showed that the minimum distance 
of the quantum code is determined by its enumerators. In a strong 

connection between a big class of quantum error-correction codes, which can 
be seen as an analog of classical linear codes, and self-orthogonal codes over 
GF(4) was found. 

A quantum code, say Q, of length n and dimension denoted as 
{{n,K)) code, is a dimensional subspace of the Hilbert space C^". Dur- 
ing transmission through a channel a code word can be altered by an error. 
In general a quantum channel error is an operator E acting on C^". If a 
codeword v G Q has been effected by an error E then v becomes E-v. 

Here we consider one of the most popular channel models - the com- 
pletely depolarized channel. In this model any error operator E can be 
represented in the form of the tensor product of Pauli matrices and two by 
two identity matrices 



E = CJi® ■ ■ ■ ® CTn, 

and fjj G {±l2,±CF^,±az,±CFxCFz] , where 



1 

1 



and CTz 



1 
-1 



The error cTx is called the flip error, and the error az is called the phase 
error. Error operators acting on C'^" form the extraspecial group of order 
22n+i 'pj^g number of nonidentity matrices in the tensor product (|l|) is 
called the weight of E and is denoted by wt{E). A quantum code can detect 



an error E if and only if |1C] 



V,' Em, 



0, 



where Vj and Vj are any two orthogonal vectors from the code. A quantum 
code has the minimum distance w and is denoted as {{n, K,w)) code if it 
can detect any error of weight less than or equal to w — 1. 

Linear codes play a special role in the classical coding theory. Similarly in 
the quantum coding theory there exist the, so called, stabilizer codes. They 
can be considered as an analog of the classical linear codes. The formal 
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definition of the quantum stabilizer codes is following. A quantum code Q 
is called stabilizer code if there exists a subgroup £ of the extraspecial group 
such that for E^^ = v for any v G Q and any E £ £. In other words, Q forms 
an eigenspace of £. If the group £ has the order 2* then dim{Q) = 2"^* and 
so the dimension of a stabilizer code is always equal to a power of 2. A nice 
property of the quantum stabilizer codes is that they are strongly related 
to classical self-orthogonal codes over G-F(4) ||2j[^. Namely, there exists a 
quantum stabilizer [[n, A;]] = {{n,2^)) code Q with the minimum distance 
w if and only if there exists a group self-orthogonal code C of length n, 
cardinality |C| = 2""^^ such that w = min{wt(v) : v G C"*" \ C}, where C"*" 
is the dual to C code with respect to the trace inner product. The trace 
inner product, denoted as *, of vectors v and u is defined as follows 



V * u = tr j ^ ViUi 



where the bar denotes conjugation in GF(4). If code C, assosiated with 
quantum stabilizer code, is linear over GF{4), then Q is called linear stabi- 
lizer code, or just linear quantum code. 

A Gilbert- Varshamov type bound was obtained in [^. 

Theorem 1 ^ There exist quantum codes of length n and dimension K 
such that 

loqK w ^ / 
^ >1 logs- H ' 



n n \n 

where H(x) = —xlogx — (1 — x) log(l — x) is the binary entropy function 
and log is the base 2 logarithm. 

The goal of the present paper is to obtain upper asymptotic bounds on 
the minimum distance of an arbitrary quantum code using linear program- 
ming approach and some strengthenings of the bounds for linear quantum 



codes. For earlier known upper bounds see Laflamme and Knill [10|, Rains 
[16|, Cleve B. 



The paper is organized as follows. In section 2 we prove the key theo- 
rem that allows us to apply the linear programming approach to obtaining 
bounds for quantum codes. In section 3 we give a short description of prop- 
erties of Krawtchouk polynomials that we will need later. In section 4.1 and 
4.2 we use the key theorem to obtain Singleton and Hamming type bounds 
for arbitrary quantum codes. In section 4.3 we strengthen the Hamming 
type bound for the linear quantum codes. 
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2 The key inequality 



In this section we obtain the key inequality allowing us to reduce the problem 
of upper bounding the size of codes to a problem of finding polynomials with 
special properties. 

Like in the case of classical codes one can introduce the notion of enu- 
merators of a quantum ((n, K)) code. An ((n, K)) quantum code Q has two 
enumerators |17|,|14| 

1 " 

wt{E)=i 

1 " 
wt{E)=i 

where E is an error operator in C^". It is shown in ||l^,||l^ that < Bi < 
B^- , Bq = Bq = 1 and the values Bi and B^- are connected by Mac Williams 
identities, Bi = ^ ^]Lq where Pi{t) is the 4— ry Krawtchouk i-th 

polynomial and S = j=o ^j' ■ 0] ' [0 ^^^'^ shown that the minimum 
distance of Q equals the maximum integer w such that Bi = B^,i < w — 1. 

It is easy to check that K = — — — = and thus we are interested in 
estimating the value w for given values n and S. 
Let f{x) be a polynomial of degree at most n, 

n 

f{x) = Y,m{x). 

i=0 

Let, moreover, all the coefficients fi be nonnegative, /(x) > for x = 
0, ... ,11! — 1, and f{x) < for x = w, . . . ,n. Then 

to— 1 n 

Sj2fiBi<SY.fiBi 

i=0 i=0 
n -, n n n 

= c/^ E Efnu) = E E f^j) 

i=0 j=0 j=0 i=0 

n w—1 to— 1 

= E/0')^/ < E f(3)B^ = E fim- 

j=0 j=0 j=0 
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Thus, 

fU) 

< max — - — . 

j=0,...,w^l fj 

We formulate now the result as a theorem. 
Theorem 2 Let Q be an {{n, K,w)) quantum code. Let 

n 

fix) = mix) 

1=0 

be a polynomial, fi > 0, and fix) > for x = 0, . . . ,w — 1, and fix) < 
for x = w, . . . ,n. Then 

r.^ S 1 fij) 

K < — < — max . 

- 2" - 2" j=o,...,w-i fj 

3 Krawtchouk polynomials 

Here we survey some properties of Krawtchouk polynomials. We consider 
the quaternary case. In this case the Krawtchouk polynomials are defined 
as follows: 

p.(.)^i:(-iF3-(;) (•;:;). 

Every polynomial of degree at most n has a unique expansion in the 
basis of Krawtchouk polynomials. If a polynomial fix) has the expansion 

t 

fix) = Y,f\P,ix), 

1=0 

then 

n 
j=0 

The following property (see [^]Chap. 5, Exercise 41) is important: 
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The Krawtchouk polynomials satisfy a recurrent relation, 

{i + l)Pi+i{x) = (3n -2i- 4x)P,{x) - 3{n - i + l)P^-l{x). 
Some useful values of the Krawtchouk polynomials: 

Po{x) = l, Pi(x) = 3n-4x, 2^2(2;) = 16x^ -8x(3n- 1) +9n(n- 1) 

Pii0) = 3' 



n ' 



We will also need Christoffel-Darboux formula for binary Krawtchouk 
polynomials 

Pt+i{x)Pt{a) - Pt{x)Pt+i{a) (3) 
2(a — x) fn^ ^ Pi{x)Pi(a) 



t + l \t 



i=0 



4 Upper bounds 

In what follows we present several bounds derived from Theorem |2|. 

4.1 Singleton type bound 

Choose in Theorem 



n { n—x 

X, 



/(x)=4--+in(i--)=4 



n—w+l \n-w+V 



Then 



j=0 



n ( n-j \ 

j=0 V«i-1/ 

/n— a;\ 
_ \w-V 

C-l) 
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Now, 



rf'T^ — ■^^ — An~w+1 (n-w+l) 
' \-^) — f — ^ (n-x\ ■ 

\w-ll 

Considering the ratio 

/ s / n—x \ ln—x—l\ 

f'yX) _ \n-w+l) ^ V w-1 I 
X 



r{x + l) r-^) " C^-^-l) 

n — X — w + 1 
w — X — 1 ' 

we find that it is greater than 1 \i w < {n + 2)/2. So, the values r[x) are 
decreasing, and we have the following theorem. 



Theorem 3 



1 /(O) r,n-2w)+2 



K < — --r^ = 2 

- 2" fo 



Note that though this bound has been already derived in |10|,[16|, the 
proof presented here is different and so could be of interest. 

4.2 Hamming type bound 

Let e = iw- l)/2. Define ^ = {Pe{x)f. 
Lemma 1 

Pi{x)Pj{x) = 

t 1:7, ^ Un - A /2A: + 2. - z - j\ 

to hVk + 2s-^-J)[ s )[ k + s-j ) 

Proof The proof of the lemma is a straightforward generalization of the 
proof of the similar expression in the binary case (see e.g. |jll|(A.19)). □ 
Using the lemma, we get 

This yields 

m = 
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n n n—k 

EEE 

j=0 k=0 s=0 



k 

2k + 2s-2e, 



n n—k 

EE 

fe=0 s=0 



k \ (n-k\ (2k + 2s-2e\ 2e-2s-kos 



2k + 2s- 2e, 



k + s — e 



s=0 



e-x/2 

= 4" E 

s=niax{0,e—x} 

Taking into account that 



X 



2x + 2s-2ej 



. s I \ X + s — e I 



n 



H - +0 - , 



n 



and denoting ^ = x/n, u = s/n, and r = e/n, we get 

X \ In — x\ (2x + 2s — 2e ' 



-log 
n 



2x + 2s- 2e 

2^ + 2v- 2t 



k + s — e 



22e— 2s— fcgs 



+ (1-0^ 



1-^ 



+ z/ log 3 + O 



Taking derivative of the last expression, we get 



d (CH (2£±f^) + (1 - OH (^) + . log 3 + O (i)) 



21og 1 



2i + 2v- 2t 



+ log 1 



V 



2 log 



1-e 



log 



1-c 



2^ + 2z/ - 2t 
+ log3. 



It is not difficult to check that this function has only one root on the interval 
max{0, T — <v <T — ^/2. Let a(T, ^) be the root of the function in the 
interval. Then we can estimate Mog/(x) as follows 



1 



n 



=2+^+^H 



logf{x) 
2e + 2a(r,e)-2r 



(4) 
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+(1 - OH i^}^) + air, log 3 + O (- 



To obtain an estimate on fx we need bounds on values of Krawtchouk 
polynomials. We will also need a bound on r^, the smallest root of Pe{j)- 
For t growing linearly in n and e = rn (see e.g. P]) 



= ^ = ^ - - l^3r{l-r) + (5) 

The next observation is a generalization of a similar fact in the binary 
case iPl. 



Lemma 2 For C < 



-logPe(x) = -log (")e=^+ (6) 



n n 



Sog f ^ - - + ^'^ - - ■'^'^ - '^-'^ - -I) + 0(1; 

V 6(1 -z) In 



Proof Like in the binary case (see, e.g. [12, Lemma 36]) one can show that 
for C <Ct 

PtU + 1) ^ 3n - 2j - 4t + V(3 - 2j - Aty - 12j{l - j) ^ 1 
Ptii) 6(1 -j) n 



Taking logarithm on both sides, applying this recursively to i-t(O) 

t)' 



'^")3*, and approximating the sum by the integral, we get the claim. 



□ 

Notice, that the integral in Lemma || can be expressed explicitly, namely 
, / 3 - 2z - 4r + Ji?, -2z- 4r)^ - 12z(l - z)\ , 

[ e^fe ) 

= —z log 6 + (1 — z) log(l — z) + z log(c + t) 

log(3 + a-^z-2t) + log(a - 2) - log 2 

log(6 + 2a-a^ - lOz + 2az - (a - 2)t), 
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where a = 3 — 4r, c = a — 2z, t = \/ 6^ — Vlz(\ — z). Hence 

-log/, = 2rlog3 + 2/7(r) (7) 

n 

+ (-Z log 6 + (1 - z) log(l - z) + z log(c + t) 

log(3 + a - 8z - 2t) + log(a - 2) - log 2 



log(6 + 2a - - lOz + 2az - (a - 2)f 



n 



Using Lemma |2|, we can estimate ^ log as follows 

-log/, = 2rlog3 + 2/7(r) 

n 

+ f log { 3-2^-''r+V(3-2.-4r)-^-12.(l-.)\ _^ ^ 1 

JO \ 6 — 62; j n 

Using (0) and (Q), we get a Hamming type bound 

logJC 
n 

2e + 2a(r,0 -2r 



< max {l + £ + £^H ^ 
~ 0<C<5 IV i 

+(1 - 0^ ( yi^) + «(^' 6 log 3 - 2t log 3 - 2H{t) 

' log [^3-2.-4r + V(3-^2.-4r)^-12.(l-.)^ 

where 5 = 2t = w/n. Analytical computations with Maple show that this 
function achieves its maximum at ^ = for any 5 < £,e- From this condition 
and ^ it follows that (5 < in the interval < 5 < a, a ~ 0.34. 

Note that when 5 > the coefficient /, can be equal to zero and the 
bound tends to infinity. 

So we conclude that the conventional Hamming bound is valid when 
5 < 0.34. 

Theorem 4 If Q is an {{n,K)) quantum code and (5 < ^g; then 

< 1 log 3 — H 



n ~ 2 V2 
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Notice that would the Hamming type bound be proved for the total 
range, the value ^ would be equal approximately to 0.38 when log(i^)/n = 0. 
This is worse than the bound 5/n <^ derived by Rains in 

We labeled quantum Hamming bound on figure 1 by "H" . On the figure 
1 we also present so called the first and the secon bounds (bounds LPl and 
LP2), though for the time being it is not proved that they are valid. We 
shall prove in the next section that in the binary case the conventional first 
linear programming bound is valid for quantum codes. 



4.3 The First Linear Programming Bound 

In the present version of the paper we confine ourselves by the binary case. 
Work on the quaternary case is under process. 

We formulate the problem the same as in the quaternary case. We have 
numbers Bi and B^- such that < i?j < B^^, Bi and B^- are connected by 
binary MacWilhams identities, and (E"=o^/") (Er=o^i) = Let w be 
the first integer such that B^-i < B^_i. Denote S = Y^^=o Like in the 
quaternary case we are interested in an upper bound for the value S given 
n and w. 

Let 

fix) = {Pt+iix)Ptia) - Ptix)Pt+iia)}^ . 

a — X 

This polynomial allows to get the so called first linear programming bound 
for classical codes Denote by x^*^ the first root of the binary polynomial 
Pt{x). To get the first linear programming bound one has to choose - = 

i - ^|W^) + 0(1) and .r^^ <a<x?,T^ = -I. 
Using Christoffel-Darboux formula (^), one can write 

' Pi{x)P,{a) 

i=0 

'- {Pt+i{x)Pi{x) + Pt{x)Pi{x)} . 
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Using identity 



P,(x)P,(x) = ^P,(x) 



r+s-j I I r-s+j 




j=0 \ 2 / V 2 

we get 

+EfS)(iJ±.l^.(-) 

j=Q\ 2 / \ 2 



+ 



^- J W J 

t+i-j t~i+j 



Suppose now that t and j are even numbers (in other cases arguments 
are similar). Estimate the sum in the last expression by the term for i = t. 
Then we have 

2p2(a) (n-j\ (j\ 



fj> 



t + 1 



Hence 



fix) ^ ^ {Pt+i{x)Pt{a) - Pt{x)Pt+ii a)} 
t + l {Pt+i{x) + Pt{x)y 



2 



2{a-x) 

2 2 

Replacing Pt{x) by P.(i)|i| and Pt+i{x) by P,(i + 1)^, we get 

f{x)^ t + l {{-,)PAt + i) + OPAt)r 



U ^ 2(a-x) (r^)(|)Q' 



2 

a;-' 
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Now using the relation 

Pxit + l) _ f^_^Qf'^\\ n-2x + - 2x)^ - 4x(n - x) 



Px{t) \ \nj J 2n-2x 

and replacing Px{t) by Pt{x)j^, we get 



fn\'^lD f+\2 j {n~t){n-2x+^/ {n-2x)'^ -ix{n-x) 
f{x) t + 1 I (t+l)(2n-2x) +^ 



< 



" 2{a-x) (r^-)(s)n' 

2 2 



{n-t){n-2x+^y(n-2x)'^-4x{n-x))+{t+l){2n-2x) 



1 -ft(a;)^ ^ 2(n-x) 



2{a-x) + 

2 2 

Using the estimate 

1 log P,(.) = + f log ( i-^^wa-2rf-4.a zi)) + 

n JO \ 2 — 2z / 



where ^ = x/n and r = t/n, we have 



li fix) ... ^ /"^ / 1 - 2r + - 2t)^ - 4z(l - z) \ , 
hm - log = 2i?(r) + 2 / log — -^^ ^ dz 



n^oo n fx Jo \ 2 — 2z 

Analytical computations with Maple show that for r > 0.11 this function 
achieves its maximum at ^ = 0. Since ^ = ^ = ^ ~ V'^i^ ~ ''") have the 
following theorem. 

Theorem 5 If 6 < 0.1865 then 
1 



-logS<Hi^--^6{l-6) 

Note that we are interested in an estimate of w for all S such that - log S > 
0.5. Note also that when 6 < 0.1865 then -log 5 ^ 0.501. So for ah 
1 > ^logS > 0.501 the conventional first linear programming bound is 
valid. 
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4.4 Strengthening in Linear Case 

Let now Q be a stabilizer quantum code with associated self-orthogonal code 
C over GF{4). Recall that if dim{Q) = 2^ then |C| = 2""'=. 
A generator matrix of C can be written in the form 



G 



4i ^3 



(8) 



where Bj is a binary matrix, Aj is an arbitrary matrix, and uj is an element 
of GF(4) of order 3. We will say that G defines a code of type 4''o2''i. If 
in some / coordinates a group code, say D, of length n contains only 0-s 
and a-s, a G GF(4),a 7^ 0, then we will say that D is a mixed code of 
lengths / and n — l. Note that if C is a code of type ^^^2^^ of length n then 
corresponding quantum stabilizer code Q is an [[n, k'\\ = [[n, 2n — 2kQ — fci]] 
code. 

Lemma 3 The minimum distance w of Q is not greater than the minimum 
distance of the optimal group mixed code of lengths ki and n — k^ — ki and 
cardinality 22"-4A;o-2fci _ In particular 

i) if ki = then w is not greater than the minimum distance of the 
optimal group code of length [n + k)/2 and cardinality 4'^. 

ii) if ki < ^"-""^^o = 2k then w is not greater than the minimum distance 
of the optimal group code of length and cardinality 2'^^~^^ . 

Proof 

Let Gc be a generator matrix of the code C written in the form (^). 
Since G C C"*- we can append some rows to the matrix Gc to get a generator 
matrix of C^. That is 



Gn-L 



G' 

Let us call the code generated by the matrix G' a complementary code of C 
to C^. It is clear that the minimum distance of a complementary code has 
to be not less than w = min{wt(v) : v G C"*- \ C}. Due to the structure (^) 
of the matrix G we can make elements of G' on the first /cq positions be equal 
to and elements on the next ki positions be equal to or lij. So G' will 
have the form[0 Di D2] where Di is an 2n — 4/i;o — 2ki x ki matrix consisting 
from 0-s and uj-s and D2 is is an arbitrary 2n — Ako — 2ki x n — ko — ki 
matrix. 
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i) Follows from the previous. 

ii) In the case ki < ^^^^ the matrix [0 Di D2] can be transformed to 
the form 

^0 Ai Bi 
^ ^2 

where Ai is an ki x ki matrix consisting form 0-s and lo-s, Bi and B2 are 
arbitrary ki x n — ko — ki and 2n — 4/co — 3ki x n — kQ — ki matrices. Since 
the subcode with the generator matrix [0 B2] has length n — ko — ki, 
dimension 2n — 4/co — S/ci and its minimum distance has to be not less than 
w the assertion follows. □ 
Let D he a mixed code of lengths / and n — I and dimension k. To get 
bounds for Q we have to get bound for the classical mixed code D. Plotkin 
and Hamming type bounds for mixed codes can formulated as follows. 



Lemma 4 



d < 



1 



2>' - 1 



{12 



3(n - /)2*^-2). 



Proof 

Let M be an array of all codewords of D. The total number of nonzero 
entries of M is equal to / • 2*^"^ + (n — /) • 3 • 2^~'^. Since the number of 
nonzero codewords in D is 2*^ — 1 we get the assertion. □ 



Lemma 5 



were e 



d-l 



EE( 

1=0 j=o vJ 



< 2 



2n-l-k 



Proof 

Let F be the Cartesian product of GF(2)' and 0^(4)""'. The volume 
Ve of a sphere of radius e in F is 




1=0 j=0 



I An~l 



The number of codewords of D can not 



d-i 



and assertion follows. 



□ 



and the volume of F is 2 4 
exceed the value V/Ve where e 

Combining Lemma ^ and Lemma ^, we get a Plotkin type bound for a 
code of given type 4'^'' 2*^1. 
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Theorem 6 



, n + A;3-4*=-i ki 

d < ; + —-, . 

- 2 4^-1 2 4^ - 1 



Combining Lemma |5| and Lemma ^, we get Hamming type bound for a 
code of given type 4'^''2''=i. 



Theorem 7 



n — k 

4— 



were e 



i=o j=0 
d-l 



3 



Let now ki = 0. For example /ci = if the code C is linear. In this case 
a complementary code is not a mixed code. So applying the asymptotic 
version of the second linear programming bound to this code we get a 
bound on the minimum distance of the code Q. One can see (bound S 
on fig.l) that this bound is better than the first linear programming and 
Hamming bounds on some interval. 

Let ki > and n — ki < k < ki/2. Then according to Lemma ^ we can 
estimate the minimum distance of Q as the minimum distance of a group 
of length "-+^"^1 and cardinality 2'^^~^'^ . For example, using the asymptotic 
bound from we get bounds for ki = 0,ki = 0.2, fci = 0.5 on figure 2. 
So for small values of ki we still have some strengthening of the first linear 
programming and Hamming bounds on some interval. 
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